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Abstract 

We consider the quantum field theory of a single, immobile, spin S hole coupled to a two- 
dimensional antiferromagnet at a bulk quantum critical point between phases with and without 
magnetic long-range order. We present an alternative derivation of its two-loop beta function; the 
results agree completely with earlier work (M. Vojta et al, Phys. Rev. B 61, 15152 (2000)), and 
also determine a new anomalous dimension of the hole creation operator. 
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Recent papers JT|, have introduced the following model Hamiltonian for a single non- 
magnetic (Zn or Li) impurity in a two-dimensional ci-wave superconductor or spin-gap insu- 
lator (see pi for a review and experimental motivation): 



Tt = - ^foS a (f) a (x = 0) 
Us = I d d x 
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We have written the Hamiltonian in d spatial dimensions, and S a (a = 1, 2, 3) are spin S 
operators of a magnetic moment that is postulated to be present near the impurity (the case 
of physical interest has S = 1/2); these operators obey the SU(2) commutation relations 

[S a , Sp] = ie a/ 3^S y (2) 

and S a S a = S(S+ 1). The field 4> a (x, t) represents the local orientation of the antiferromag- 
netic order parameter at spatial position x and time t; its canonically conjugate momentum 
is 7r a (x,t), and hence 

[<f> a (x, t), Trp(x', t)] = t6 af3 5 d (x - x') (3) 

This theory has a bulk quantum critical point at s = s c between a phase with magnetic 
order (s < s c , (<p a ) ^ 0), and a symmetric phase with a spin gap (s > s c , (<p a ) = 0). 
We are interested in the spin correlations of 7i for s close to s c , and in the vicinity of 
the impurity at x — 0. As discussed in JT], 0, universal aspects of these correlations are 
associated with a renormalized continuum theory of Ti defined in an expansion in e = 3 — d. 
This renormalization involves the familiar bulk renormalizations which are insensitive to the 
impurity degree of freedom 

<p a = \fz$ Ra ; go = i?<f 4 9 (4) 
and new 'boundary' renormalizations associated with the impurity spin 

s« = Vz i s Ita ■ 7o= r r ( 5 ) 

\IZZ'S d+1 

Here [i is a renormalization momentum scale (we set the velocity c = 1), = 
2/\r{d/2)(ATi) d % and S d = T(d/2 - l)/[47r d / 2 ]. The renormalization constants Z, Z A were 
computed long ago H; their values in the minimal subtraction scheme to order g 2 are 

Z = 1 ~lUe ! Z * = l + ^ + {-^-We) 9 ■ (6) 



The boundary renormalizations were computed to the same order in [I]. [2|] : 

7=1 1 , Z 7 = H 7 g (7) 

e e be 

This paper will rederive the above results by a new method which also yields a renormal- 
ization constant for the hole creation operator. Furthermore, the present approach, unlike 
that of 0, has the advantage of being formulated entirely in terms of perturbation expan- 
sion which has a Wick theorem, and can thus be presented in conventional time-ordered 
Feynman diagrams. 

We will identify the spin S a with that of a hole, with creation operator t/>t, that has been 
injected into the antiferromagnet. So 

S a = ^ a L^ b (8) 

where a, b take the 25 + 1 values — 5, . . . S, and the L a are the (25 + 1) x (25 + 1) angular 
momentum matrices associated with the spin 5 representation. The hole operators obey the 
anticommutation relation 

^ b + Ml = S ab (9) 
So the remainder of this paper will consider the Hamiltonian 

= X^J a + H (f> - ^lL^lk<f>a{x = 0) (10) 

We will only look at the Hilbert space with a single hole, and A, the energy of this hole is 
an arbitrary positive number. 

We now consider the renormalization of TC^. The standard procedure suggests the pa- 
rameterization 



ifj a = JZ h il) Ra ; 70 = —. ~2 7- (11) 



Zh\ ZSd+i 

It is important to note that despite the relation (|]), the renormalization of the spin S a is 
not the square of the renormalization of ip a , Z' 7^ Z\\ bringing the two Fermi operators to 
the same spacetime point introduces a composite operator renormalization which invalidates 
such a relation. Instead, the relationship between the two renormalization schemes emerges 
by comparing the renormalization of 70 in and (|TTD; consistency of these relations de- 
mands 

Z\Z* = Z*Z' (12) 



We will now compute Zh and Z 7 by completely standard field theoretic methods, and verify 
that their values and (0) satisfy ([12]). 

The Feynman diagrams for the renormalization of two-point ip Green's function are shown 
in Fig [I]. As an explicit example, we display the computation of the simplest one-loop graph 
in Fig [l|a: 

rt \ 2 C / C . t ddk rdu' 1 1 
(la) = 7 5(S + 1) 



s r 

l 2 oS(S+ l)fJ Q 



(2n) d J 2ir (u/ 2 + k 2 ) {-i{oJ + u') + A) 
S H r°° k d ~ 2 dk 



-iuj + k + A) 



A^-iu + X^SiS + l) 



~ + N/2 + 0(e) 



(13) 



where = fi € (—iu) + \)~~ e Z 2 /(Z 2 Z). In the last step, the integral was evaluated in dimen- 
sional regularization. The constant K = —0.8455686701969 ... is a consequence of phase 
space factors and will eventually cancel out of our final results. The remaining diagrams can 
be evaluated in a very similar manner: the frequency integrals are performed first, followed 
by integrals over the radial momenta. The results for the two-loop diagrams in Fig [I] are 



(16) = Ai(-iu + X)^S 2 (S + l) 



1 1 - 
2? + ^2~ e 



+ O(e ) 



(lc) = ^(-^ + A) 7 4 5(S+l)(5 2 + 5-l) 



1 



-1 + 2N 



+ O(e ) 



(14) 
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Turning to the renormalization of the vertex 70, the Feynman diagrams are shown in 
Fig 0. Evaluating these as above we obtain 

1 



(2a) 
(26) 
(2c) 
(2d) 

(2e) 
(2/) 



7o^ 7 2 (5 2 + 5-l 
7o^7 4 (^ + 5-l) 2 



1 - K/2 + 0(e) 



3 + K 
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2e 



+ O(e ) 



7o A 2 7 4 (S-l)(S + 2)(S 2 + S-l) 

1 2 + K 



1 



7oA 2 7 4 (5 2 + 5-l) 2 
7oA 2 7 4 ^+l)(5 2 + 5-l) 
-7o 



2e 



2e 
+ O(e ) 



O(e ) 
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7 0(5 + 5 - 1/3) 



Z 2 Z 



(15) 



The two-loop expression for the boundary renormalization constants follows immediately 
from the results (|I3~1 ,14>15). Demanding cancellation of poles in e in the expressions for the 



renormalized vertex and ip Green's function at external frequency — iu + A = /x we obtain 
Z h 



(S - 1)S{S + l){ S + 2) + S(S + 1) 



2e 2 



2e 



(S 2 + S- 3)(S 2 + S-1) (S 2 + S - 1) 



2e 2 



2e 



+ 97 



,tt 2 (S 2 + S- 1/3) 
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(16) 



It can be checked that ([ID and (0) satisfy ([12]). 

The validity of (|12|) implies that the beta function for the coupling 7 is the same as that 
in §. Using either or we obtain 



■? + 7 3 -7 5 + ^ + ^ 2 + ^-l/3). 



(17) 



The anomalous dimension of the ip a field at the quantum critical point also follows from ( |T6"P 



Vh = Pin) 



drj 



5(5+l)( 7 2 -7 4 ), 



(18) 



while, as in 0, the anomalous dimension of the spin field, S a , follows from (^|,[7[): 

V = 2( 7 2 -7 4 ). 

For completeness we also note the beta function for the coupling g which follows from 

ll£ 2 23<7 3 



-eg 



(19) 



(20) 
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The stable fixed point of the beta functions ( |TT| , |2"0"D has g 7^ and 7^00. Evaluating 
' [18|) at the fixed point of the beta functions 0, we obtain 

\ ( 5 tt 2 {S 2 + S- 1/3)' 



Vh = S(S+l) 



t 2 + 0{t A 



(21) 
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(jf = 2r]f l /[S(S + 1)] at this order). This anomalous dimension implies that the Green's 
function G = (ip a ipl) obeys 

G(u) ~ (A - Lu)~ 1+rih . (22) 

The equations QlBU lS l pTD are the main new results of this paper. Unfortunately, the order 
e 2 corrections in ([H]) are rather large: this suggests that truncating the asymptotic series 
for rjh at order e probably gives the most reasonable estimate for its numerical value. 

There is also an unstable fixed point at which the bulk interactions vanish (g = 0). As 
shown in 0, 7/ = e exactly at this fixed point, and here we find that rj h = S(S+l)e/2+0(e 3 ). 



There appears to be no general reason for the higher order terms in r/h to vanish. The g = 
fixed point can also be studied in a large N theory j2], and the N = oo results are rj' = 1 
and i] h = 1/2. 

The physical motivations and implications of the above results are discussed in a separate 
paper [||]: there we argue that the anomalous dimension r/h characterizes photoemission 
spectra of mobile holes in two-dimensional antiferromagnets and superconductors in the 
vicinity of points in the Brillouin zone where their dispersion spectra are quadratic (i.e. 
near energy minima, maxima, and van Hove singularities). The intensively studied (vr, 0), 
(0, 7r) points (the anti-nodal points) in the high temperature superconductors are prominent 
examples. 
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FIG. 1: Diagrams contributing to the ip fermion self energy. The full line is the fermion propagator, 
while the dashed line is the (j) a propagator. 
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FIG. 2: Diagrams contributing to the renormalization of the coupling 7. The full circle is the 
interaction g 
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